function getCellCentreGrid(self)

%Nodal Grids
nodeX = zeros(self.ncx+1,1);
for ii = 1:self.ncx
    nodeX(ii+1) = nodeX(ii) + self.hx(ii);
end

nodeY = zeros(self.ncy+1,1);
for ii = 1:self.ncy
    nodeY(ii+1) = nodeY(ii) + self.hy(ii);
end

nodeZ = zeros(self.ncz+1,1);
for ii = 1:self.ncz
    nodeZ(ii+1) = nodeZ(ii) + self.hz(ii);
end

%Cell centre grids
ccX = nodeX(1:end-1) + self.hx/2;
ccY = nodeY(1:end-1) + self.hy/2;
ccZ = nodeZ(1:end-1) + self.hz/2;


[CCx,CCy,CCz] = ndgrid(ccX,ccY,ccZ);
self.CC = [CCx(:),CCy(:),CCz(:)];


